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ABSTRACT 


Some  battlefield  models  have  a  component  in  them  which  models  the  time  it  takes 
for  an  observer  to  detect  a  target.  Different  observers  may  have  different  mean  detection 
times  due  to  various  factors  such  as  the  type  of  sensor  used,  environmental  conditions, 
fatigue  of  the  observer,  etc.  Two  parametric  models  for  the  distribution  of  time  to  target 
detection  are  considered  which  can  incorporate  these  factors.  Maximum  likelihood  esti¬ 
mation  procedures  for  the  parameters  are  described.  Results  of  simulation  experiments 
to  study  the  small  sample  behavior  of  the  estimators  are  presented. 
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I.  INTRODUCTION 


Some  battlefield  models  have  a  component  in  them  which  models  the  time  it  takes 
for  an  observer  to  detect  a  target.  A  common  model  is  that  the  time  to  detect  a  target 
is  a  random  variable  having  a  distribution  /^(l  —  [Ref.  !]■  However  different  ob¬ 
servers  may  have  different  mean  detection  times  due  to  various  factors  such  as  the  type 
of  sensor  used,  environmental  condition,  fatigue  of  the  observer,  etc.  In  this  thesis  we 
will  consider  two  pa-ametric  models  for  the  distribution  of  the  time  to  detection  which 
can  incorporate  these  factors. 

Suppose  there  are  M  observers,  such  as  tank  crews.  Observer  i  is  presented  with  A' 
targets.  In  Chapter  2.  observer  /  has  a  parameter  A,  which  reflects  the  ability  of  the 

observer  to  detect  a  target.  The  parameters  of  the  M  observers.  T,.  A2,  A} . -tM  .  are 

assumed  independent  having  a  common  gamma  distribution.  Given  A,,  the  detection 
times  of  observer  /  are  conditionally  independent,  having  Weibull  distributions.  In 
chapter  2  the  parameters  of  the  conditional  Weibull  distributions  are  assumed  known 
and  interest  is  in  estimating  the  parameters  of  the  gamma  distribution. 

In  chapter  3  observer  /  has  explanatory  variables  x.,,  x..,,  x,3 . xrp  relating  to  his 

/'■  target  representing  factors  which  influence  his  time  to  detection.  The  detection  times 
for  the  observers  are  independent  random  variables  having  Weibull  distributions.  The 
location  parameter  of  the  Weibull  distribution  for  the  detection  time  of  the  j:h  target  by 
the  C  observer  is  of  the  form  ur  =  c-  -  where  x  />  =  ^x  ,/?*  ;  the  shape  parameter  is  of 

l 

the  form  cC  Interest  is  in  estintation  of  {/C }  and  {T(- 

In  both  chapters  the  data  for  the  observers  can  be  censored.  The  i"  observer  has  an 
opportunity  time  0,  to  detect  his  fh  target.  If  the/"  target  is  not  detected  in  time  O  .  the 
time  to  detection  is  censored.  Data  for  the  i'h  observer  consist  off  possibly  censored  ) 
times  to  detect  the  X,  targets. 

Chapter  2  and  3  present  iterative  procedures  based  on  maximum  likelihood  to  esti¬ 
mate  the  parameters.  Results  from  simulation  studies  of  small  sample  size  behavior  of 
the  estimators  are  given. 

These  models  and  estimation  procedures  should  be  of  eventual  use  in  the  Army 
MANPRINT  program,  an  objective  of  which  is  to  better  understand  the  human  con¬ 
tribution  to  battlefield  performance  [Ref.  2|. 
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II.  A  HIERARCHICAL  MODEL  FOR  TIMES  TO  DETECTION 


A.  MODEL 

Suppose  there  are  M  observers  who  attempt  to  detect  targets  on  a  battlefield.  Ob¬ 
server  i  has  a  parameter  A,  ,  i=  1.2,3,....M.  which  reflects  the  ability  of  the  observer  to 
detect  a  target.  Observer  i  is  presented  with  .V,  targets.  The  targets  are  presented  one  at 
a  time.  Let  U,.  be  the  time  it  takes  for  observer  i  to  detect  target  j.  Assume,  given 
A,  =  8  .  the  l\,  j  =  1,2, 3,. ...A',,  are  conditionally  independent  random  variables  with 
Weibull  distributions 

Pil'ij  <  1 1  A,  =  8}  =  1  -  exp  {  -  Odin,/' *'}  t  >  0  (2.1) 


independent  of  other  observers. 

Further,  assume  the  parameters  A„  i  =  1.2.3 . M.  are  independent  identically  dis¬ 

tributed  having  a  Gamma  distribution  with  density  function  as  follows: 


gtf)  = 


yly.8Y 

ri;) 


-It3 

e 


(2.2) 


where  8>0.  The  A  variations  are  introduced  to  represent  individual  observer  difier- 
ences.  Now  let 


Z;  =  In  C  . 

v 

It  follows  that 

P{Zij<z\Al  =  8)  =  P[i:ij<e!\Al  =  8\ 

=  1  -  exp|  -  0[  exp(r  -  In  M,y)]f  }  ■ 

Hence,  the  conditional  density  function  of  Z.  ,  given  A,  =  8  ,  is 

Jz  i  A{:  =  9  c.xpjlr  -  In  ij.)c~  1  exp[  -  0(  expjlr  -  In  ny)e~  *'])]  .  (2.4) 

\V  lien  the  /•"'  observer  is  presented  with  his  j:l-  target,  he  gets  a  length  of  time  called 
opportunity  time  O  to  detect  it.  An  observer  cither  successfully  detects  the  target  within 
this  time  or  i>>  unsuccessful.  Data  for  the  /•"  observer  consist  of  times  of  detections  for 


the  successes  and  the  lengths  of  opportunity  times  for  the  failures.  For  each  i  = 
1.2,3 . M  ,  j  =  1,2,3...., A',,  let 


and 


Yy  =  min(  In  Uy,  In  Of 


(2.5) 


1  if  Uy  <  Oy 
0  otherwise  . 


(2.6) 


The  are  the  censored  ln-detection  times  and  A,.  is  an  indicator  of  whether  or  not  the 
ln-time  to  detect  the/*  target  by  the  /'*  observer  is  censored.  Let 

y 

=  {2J) 
7=1 

be  the  number  of  targets  detected  by  observer  /. 

In  this  Chapter  we  will  assume  {/j,  }  and  {c.}  are  known  constants.  We  are  interested 

v 

in  estimating  the  parameters  rj  and  /?0  with  a  =  t”  and  —  =  using  maximum  likelihood. 
These  parameter  estimates  can  be  used  to  predict  future  times  to  detection  for  an  ob¬ 
server  given  his  past  performance.  In  the  next  section  we  describe  a  Newton  -  Raphson 
procedure  for  solving  the  likelihood  equations,  [  Ref.  3],  [  Ref.  4],  In  the  final  section  of 
this  chapter  we  describe  a  simulation  experiment  to  study  the  small  sample  properties 
of  the  estimators.  Results  of  the  simulation  are  also  presented. 

B.  ESTIMATION 

I.  The  Likelihood  Equation  and  Maximum  Likelihood  Estimates 

Given  =  6  ,  it  follows  from  equation  (2.3)  and  (2.4)  that  the  conditional  like¬ 
lihood  function  for  observer  /  using  the  censored  In-times  y,.  is 


>•  £:  0)  =  [de^  ~ ln  V  {-]A«  exp[  -  6e^  " ln  ^ e~  . 


(2.8) 


7=1 


Let 


3 


•  I 

Si  =  Xexpfey  ” ln  f']  • 


Rewriting  equation  (2.8) 


where 


£:  6)  =  dc'  exp|^A,y[0'7  -  ln  ny)e  ’  -  fjj  exp(  -  OS,) 


=  6C-Kt  exp(  -  6Si) 


(2.10) 


Ki  =  exp-.  V AyCCv',7  -  In  nv)e  tl  -  ?,] 


The  unconditional  likelihood  for  observer  i  is 


Lt(j±.  c.  a,  y)  =  [  £i  0)  g(6)dd 


e~^dd 


(2.11) 


(2.12) 


If  C,  =  0.  then  n  (;’  +  k)  =  1.  Recall  the  parametrization  a.  =  e’’  and  —  =  .  The  un- 

*=o  _  K 

conditional  ln-likclihood  function  for  observer  /  can  be  rewritten  as 


ln  Ljiju,  L  )].  /j0)  =  ln  A)  +  —  (Q  4-  ev+  Po)  ln(S,-  +  c7) 

c,  —  i 

+  y  \n(ev+Po  +  k) 


(2.13) 


where  if  C,  =  0,  then  3)  ln(t”,'/’°  +  k)  =  0.  Since  the  observers  are  independent,  the  un- 

4-0 

conditional  ln-likelihood  for  all  M  observers  is 

v  r  C'-1  v 

ln  A  =  y  \  ln  A',-  +  en+\  ~  (Q  +  p«)  ln(S,  +  en)  +  Y  ln(^+/?0  +  m  . 

<=  1  (.  A'=0  J 


4 


(2.14) 


The  derivative  of  the  ln-likelihood  with  respect  to  rj  is 


vr 


i=  1  a-=o 


(2.15) 


The  derivative  with  respect  to  /?a  is 


r  In  /, 

CpQ 


\f 


;=) 


ln(S,  +  e  1 )  4- 


k= 0 


(2.16) 


We  are  assuming  {n,}  and  {5,}  are  known.  The  problem  is  to  find  the  maximum  likeli¬ 
hood  estimates  of  >/  and  ;  that  is.  find  rj  and  /?0  such  that  c  If1'  ^  =  0  and  ■'  ^  =  0. 


Note  that 


c  In  L 
a  i 


c  In  I. 

Wo 


\f 

V 


,,;+  p0 


i=i  L 


Cj  +  e'+h 
S<  + 


(2.17) 


Thus,  if 
equation 


-  !”  —  =  0.  then  to  solve  the  equation  =  0  ,  we  need  to  solve  the 

cfi0  ^  cr\ 
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The  derivatives  of  J[ri,  /?0)  with  respect  to  tj  and  /?0  are  as  follows  : 


(2.18) 


cj_  y  (e%  -  cy 

~  1_J  (Si  +  e")  (2.19) 

(=! 


IL 

ep0 


(2.20) 


Solving  the  equation 


c  In  L 

cP  0 


0  is  equivalent  to  solving  the  equation 


0  =  g(ri.P  0)  = 


\f 

V  ~  ln(S;  4-  ev)  + 

i=  i 


c,-i 


k= 0 


The  derivatives  of  g(>/,  /?„)  with  respect  to  rj  and  /?0  are 


(2.21) 


(2.22) 


6 


(2.23) 


M  C<  _1 


1=1  k= 0 


where  if  C,  =  0,  then  the  sum  involving  C,  -  1  is  zero.  A  Newton  procedure  to  soive  the 
cL  cL 

equations  =  0  and  -rr-  =  0  would  use  the  following  equations 

Ct]  Cp0 


o=V./?c(  + 


(n  -  t]°)  + 


Wv0W°o) 

Wo 


{Po~P°o) 


(2.24) 


0  .0.  ,  W<V°Wo)  0  cg(rj°Wl)  Ox 

0  =  g(n  ,  Po)  + - t- - (»/-»?)  + - To - (Po  ~  Po) 


c> 1 


cPo 


(2.25) 


where  tf  and  PI  are  current  values  for  t]  and  /?„.  However  note  that  if  f[i],  fi0)  =  0.  then 

cf 

~r~  =  0.  Hence,  a  Newton  procedure  that  is  more  stable  numericallv  would  use  the 

Ct] 

equations 


which  results  in 


n  7  o  r,n  Wv°W n)  o. 

0  =A>1  ■  Po)  + - - (Po  ~  Po) 


r>  ,  o  ,o,  ,  W-/*o)  ,  o, 

0  =  g(n  .  Po)  + - - (’/  -  V  ) 


Po-Po  =  - 


-Jlrj°.Pn) 

W  ,  o 
—Jn.h) 


0 

no  -v o  = 


-  /?p) 


(2.26) 


(2.27) 


(2.28) 


(2.29) 


2.  Initial  Condition 

In  this  subsection,  we  describe  a  rough  way  to  provide  initial  estimates  of  tj  and 
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fi0  to  start  the  iterative  Newton  procedure  of  the  previous  section. 

The  uncensored  random  variable  Utl  has  the  same  distribution  as 


IV::  , 

r-  r  ‘j  \e'‘ 

l‘ij  =  — ) 


where  Jl'.  is  unit  exponential  random  variable.  Let 

zij  =  In  Uy 

=  In  uy  +  e*‘  In  IVy  —e’1  In  At  ; 


then 


£[  In  IVij]  =  -  0.5772 


2 

VAR{  In  Wy]  =  — 


(2.30) 


(2.31) 


(2.32) 

(2.33) 


[  Ref.  5:  p.943  ].  An  approximation  to  the  moments  of  In  A.  is  given  by  the  first  two 

*.  i 

terms  of  a  Taylor  expansion  of  ln(.-f)  about  the  mean  £[.4,]  =  -j- 


ln(.-lf)sln(-y-)  +  — I ■— (^--s-); 

'  Or) 


(2.34) 


thus 


Let 


£[  In  A;]  s  In  y  —  In  a  ;  (2.35) 


VAR[  In  A,]2s(  y)2V’AR[AJ 

~  (  y  >  a2  (2.36) 

1 


vij  =  (zij  “ ln  m,7)c 
=  In  IVy  -  ln  . 


Hence, 
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ElVy]  =  -0.5772  -  ( In  y  -  In  a) 


(2.38) 


Ty]  =  VAR[  In  IV  A]  +  VAR\_  In  Atj]  -  2  cov{  In  IV tp  In  A,) 


=  -r+± 

6  / 


Rough  estimates  of  the  first  and  second  moments  of  V0  are  given  by 

Xf  V 
>=  i  7=1 


3/,= 


M  S, 

E2a 

<=1  7=1 


and 


M 


y  yAyi|z(/'  - in  Vi/ie  ]2 


/=’.  7=1 


v/  V 

(=1  ;*=  1 


Recall  the  parametrization 


a  =  e 


v  . 


Thus 


v  =  e  V° 


and 


£[ 1  ’;y]  =  -0.5772  —  in  y  +  In  a 

=  -0.5772 -[„  +  /?„]  + >7 

=  -0.5772  -/?,  ; 


(2.39) 


(2.40) 


(2.41) 


(2-42) 

(2.43) 


(2.44) 


(2.45) 
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(2.46) 


VAKV$-\  +  ± 


=  _E_  4.  p-  (-7+  M 
6 


Equating  the  rough  estimates  of  the  moments  of  Vv  with  the  mean  and  variance  of  VtJ 
we  obtain  initial  estimates 

i0  =  -  A/,  -  0.5772  ;  (2.47) 

V  =  -  In [A/2  -  (A/,)2  -  ]  -  k  .  (2.4S) 

'  w2  A 

If  .U;  —  (AM2  —  <  0,  we  set  r;  =  —  /?„  .  The  latter  condition  usually  occurs  in  the 

simulations  indicating  that  a  better  initial  condition  might  be  found. 
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Implementation  in  Simulation 

The  following  flow  diagram  describes  one  replication  of  a  simulation 


C.  SIMULATION  PROCEDURES  AND  RESULTS 


1.  Simulation 

All  simulations  were  carried  out  on  an  IBM  3179  G  computer  at  the  Naval 
Postgraduate  School  using  the  APL  GRAFSTAT  random  number  package  [Ref.  6]. 
Histograms  of  simulated  estimates  were  produced  by  an  experimental  APL  package 
GRAFSTAT  which  the  Naval  Postgraduate  School  is  using  under  a  test  agreement  with 
IBM  Watson  Research  center,  Yorktown.  Height,  NY.  The  simulation  is  replicated  for 
R  =  100  replications.  Each  simulation  experiment  of  100  replications  starts  with  the 
same  random  number  seed.  The  mean  bias  (M.B)  and  its  standard  error  (  S.E(M.B)  ). 
mean  square  error  (M.S.E)  and  its  standard  error  (  S.E(M.S.E)  )  from  R  =  100  repli¬ 
cations  are  computed  as  follows  : 


a 


i—  i 


R 


S.E(M.B) 


^  ((ft 

i—  i 


( 2.49a) 


( 2.49/0 


rv 


.50(0 


S.E(M.S.E)  =  /  )  {(Pi  -  P)2  -  M.S.E)2 


V 


i2.5<»/’ 


i=i 


where  /?,  is  the  point  estimate  of  the  true  value  P  for  the  /th  replication  and  R  is  the 
number  of  replications. 

The  simulation  experiment  to  study  the  sampling  properties  of  the  estimators 
of  t;  and  /?„  is  as  follows. 

1.  Give  arbitrary  constants  for  ur  and  c,  (  #i„  =  4.2  .  I,  =  0  ) 

2.  Set  the  same  opportunity  time  .  O for  M  observers,  i  =  1,2.3 . M .  /  =  1.2.3 . 

A’..  The  constant  opportunity  time  is  changed  to  give  different  censoring  levels.  The 
different  values  for  O,  are  10,  25,  40. 


12 


3.  Generate  gamma  random  variables  for  M  observers  (  A,A3,  Ay . AM  with  shape  pa¬ 

rameter  y  =  1.5,  scale  parameter  a  =  6.5  ). 


4.  Generate  independent  exponential  random  variables  with  mean  1  (  U',„  i  - 
1.2.3 M  J  =  1,2.3 . V,  ). 


5.  Compute  the  detection  time  L’u  for  observer  i  to  detect  target  j 


6.  Compare  In  l'r  with  In  0,„  then  choose  the  smaller  one  for  the  data 

}]j  =  min(  In  l  y  In  Oiy) 

and  compute  the  censoring  indicator  A,  as  in  equation  (2.6) 

7.  Compute  the  moments  .U.  and  M:  for  the  observations  that  are  not  censored  as  in 
equation  (2.40  and  equation  (2.41). 

8.  Compute  the  initial  value  for  >/  and  /?0  as  in  equation  (2.47)  and  equation  (2.48) 

6.  Lse  equation  (2.28)  and  (2.29)  to  find  new  values  for  j;  and  />„ 

10.  Iterate  the  procedure  until  the  dilfcrcnces  between  successive  values  of  >/  and  /»c 
are  small.  (  less  than  1"  "  ). 

2.  Results 

In  this  section  results  from  the  simulation  experiments  will  be  reported.  Simu¬ 
lation  experiments  were  done  for  various  numbers  of  observers  and  targets  and  values 
of  the  opportunity  time.  The  numbers  of  observers  considered  arc  M  =  5.  15.  and  3<>. 
The  numbers  ol  targets  considered  arc  .Y  =  5,  15.  3<*.  and  5c  for  each  observer.  Tor  all 
simulation  experiments,  the  u  =  4.2.  c,  =  ",  y  =  1.5  and  j  =  6.5.  I  rom  equations 
(2.42)  and  (2.44)  the  true  values  of  t]  and  /?t.  are  found  to  be  1.8" IS  and  -1.4663. 
Histograms  of  the  difference  between  the  estimates  and  the  true  values.  >/.  -  t]  and  />,.  - 
p,  are  shown  in  figures  1  through  18  of  Appendix  A.  For  Figures  1  through  6.  O  = 
1":  for  Figures  T  through  12.  O  =  25:  lor  Figures  13  through  IS.  O,  =  4i>.  The  fraction 
of  detected  level  or  uncensoring  level  (  L  C  )  for  R  replications  of  each  simulation  is  dc- 
ternuned  by  following  equation 


where  C[r)  is  numbers  of  targets  detected  by  observer  in  the  replication,  S,(r)  is  the 
number  of  targets  presented  to  observer  i  in  the  r,h  replication.  The  L'C's  for  the  simu¬ 
lations  using  the  same  opportunity  time  are  then  averaged  to  obtain  the  mean  L‘C  for 
that  opportunity  time.  For  each  Figure  of  Appendix  A  and  B,  the  mean  uncensoring 
level  for  the  opportunity  time  is  given  in  parenthesis.  Mean  bias  and  mean  square  error 
are  recorded  in  tables  1  through  9  for  each  of  the  estimators;  their  standard  errors  appear 
in  parenthesis  below.  Tables  1.  4  and  7  present  all  of  the  means  and  standard  errors  for 
different  opportunity  times.  The  other  tables  present  the  same  results  in  a  more  con¬ 
venient  fashion;  the  mean  biases  are  displayed  with  the  mean  square  error  in  parenthesis. 
For  each  Table,  the  opportunity  time(O)  and  average  uncensoring  level(UC)  for  that 
opportunity  time  are  given  at  the  top  of  table. 

•  Results  for  >/ 

The  histograms  of  >)  -  >]  are  centered  about  0  with  a  slight  amount  of 
skewness  to  the  right.  Increasing  the  opportunity  time,  which  results  in  less  cen¬ 
soring.  has  very  little  effect  on  mean  bias  and  mean  square  error.  Increasing  the 
number  of  targets  for  a  fixed  number  of  observers  has  some  tendency  to  decrease 
the  mean  square  error  and  bias.  Increasing  the  number  of  obsercers  lor  a  fixed 
number  of  targets  has  the  greatest  effect  on  decreasing  the  mean  bias  and  mean 
square  error. 

•  Results  for  />,. 

d  he  histograms  of  /»,,  -  /I  are  centered  around  0  with  some  skewness  to  the 
left.  Once  again  changing  the  opportunity  time  has  little  effect  on  mean  bias  and 
mean  square  error.  Changing  the  number  of  targets  for  a  fixed  number  of  observ¬ 
ers  also  has  little  effect.  Increasing  the  number  of  obsercers  for  a  fixed  number  of 
targets  has  the  greatest  effect  on  decreasing  the  mean  bias  and  mean  square  error. 
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Table  1.  MEAN  BIAS,  MEAN  SQUARE  ERROR  AND  STANDARD  ERROR  AT 
0=  10  AND  UC -37% 


Table  2.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0=10  AND 
I'C  =  37% 


M.B(M.S.E)  for  q 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

0.79(2.04) 

0.58(1.3) 

0.42(0.97) 

0.43(1.68) 

15  OBS 

0.19(0. SI ) 

0.19(0.72) 

0.09(0.20) 

0.14(0.23) 

30  OBS 

0.29(0.54) 

0.13(0.21) 

0.08(0.09) 

0.01(0.09) 

Table  3.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0=10  AND 
UC  =  37% 


.\/.£(3/.S.£)  for  /?0 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

-0.07(0.36) 

-0.08(0.25) 

-0.12(0.19) 

-0.05(0.14) 

15  OBS 

-O.o2i0.12) 

-0.009(0.07) 

-0.07(0.07) 

-0.02(0.06) 

3(>  OBS 

-o.03(0.05) 

-0.01(0.03) 

-0.04(0.03) 

-0.004(0.02) 
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Table  4.  MEAN  BIAS,  MEAN  SQUARE  ERROR  AND  STANDARD  ERROR  AT 
0  =  25  AND  UC=  62% 


Number  nf 

Number  nf 

h 

_ 

p 

0 

Observers 

Targets 

M.B 

M.S.E 

M.B 

M.S.E 

(S.E) 

(S.E) 

(S.E) 

(S.E) 

< 

■c mm 

2.23 

-0.06 

0.25 

mSsEm 

(0.23) 

(0.05) 

(0.07) 

15 

0.97 

-0.08 

0.20 

5 

mmmM 

(0.14) 

(0.05) 

(0.04 ) 

30 

0.47 

1.00 

-0.12 

■ 

(0.01) 

(0.14) 

(0.04) 

50 

0.47 

0.76 

-0.05 

0  13 

(0.0") 

(0.10) 

(0.04) 

(0.02) 

5 

0.20 

-0.01 

0.09 

(0.07) 

Hi; 

(0.03) 

(0.01) 

1  > 

0.15 

0.29 

-0.01 

0.07 

15 

(0.05) 

(0.05) 

(0.03) 

(0.01) 

30 

0.10 

0.20 

-0.06 

0.06 

(0.04) 

(0.03) 

(0.03) 

(0.01 ) 

50 

0.16 

■m 

-0.12 

0.06 

(0.05) 

mmwM 

(0.03) 

(0.01) 

0.17 

0.37 

-0.03 

0.04 

(0.06) 

(0.12) 

(0.02) 

(0.01) 

15 

0. 1 3 

0.18 

-0.005 

0.03 

30 

(0.04) 

(0.03) 

(0.02) 

(0.01  l 

30 

0.09 

0.09 

-0.03 

0.03 

(0.03) 

(0.01 ) 

(0.02) 

(0.01 ) 

5u 

O.o  1 

0.10 

0.02 

0.02 

(0.O3) 

(0.01) 

(0.02) 

(O.o I ) 
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Table  5.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  25  AND 
UC=  62% 


M.B(M.S.F)  for  h 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

0.83(2.20) 

0.52(0.97) 

0.47(1.00) 

0.47(0.76) 

15  OBS 

0.20(0.55) 

0.15(0.29) 

0.10(0.20) 

0.16(0.22) 

30  OBS 

0.17(0.37) 

0.13(0.18) 

0.09(0.08) 

0.008(0.10) 

Table  6.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  25  AND 
UC  =  62% 


M.B(M.S.E)  for  k 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

-0.06(0.25) 

-0.08(0.20) 

-0.12(0.18) 

-0.05(0.13) 

15  OBS 

-0.01(0.09) 

-0.007(0.07) 

-0.06(0.06) 

-0.02(0.01) 

30  OBS 

-0.03(0.04) 

-0.005(0.03) 

-0.03(0.03) 

0.02(0.02) 
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Table  7.  MEAN  BIAS,  MEAN  SQUARE  ERROR  AND  STANDARD  ERROR  AT 
0  =  40  AND  UC  =  75% 


Number  of 
Observers 

Number  of 
Targets 

h 

0 

M.B 

M.S.E 

M.B 

M.S.E 

(S.E) 

(S.E) 

(S.E) 

(S.E) 

< 

■B 

mmm 

-0.06 

0.20 

HiUH 

■■ 

(0.05) 

(0.04) 

15 

mmm i 

1.08 

-0.06 

5 

■i  - 

(0.16) 

(0.04) 

30 

0.49 

-0.12 

(0.09) 

Wam&Sm. 

(0.04) 

mmm 

50 

0.61 

1.71 

-0.05 

0.13 

(0.12) 

(0.17) 

(0.04) 

(0.02) 

S 

0.24 

0.70 

-0.01 

0.09 

(0.08) 

(0.18) 

(0.03) 

(0.01) 

15 

0.17 

0.29 

-0.01 

0.07 

15 

(0.05) 

(0.05) 

(0.03) 

(0.01) 

30 

0.12 

0.20 

-0.06 

0.06 

(0.04) 

(0.03 ) 

(0.03) 

(0.01) 

50 

0.35 

-0.03 

0.06 

(0.04) 

I 

(0.03) 

(0.01) 

5 

0.13 

0.24 

-0.02 

0.04 

(0.05) 

(0.03) 

(0.02) 

(0.01) 

15 

0.20 

SBSXi' 

-0.01 

0.04 

30 

(0.07) 

(0.02) 

(0.01) 

30 

0.16 

■P mm 

-0.04 

0.03 

(0.08) 

HUH 

(0.02) 

(0.01) 

nO 

mam 

0.09 

+  0.01 

0.02 

HffiBH 

(0.01) 

(0.02) 

(0.01) 
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Table  8.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  40  AND 
UC=  75% 


M.B(M.S.E)  for  ft 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

0.89(2.20) 

0.55(1.1) 

0.49(1.00) 

0.61(1.70) 

15  OBS 

0.24(0.70) 

0.17(0.29) 

0.12(0.20) 

0.35(2.00) 

30  OBS 

0.13(0.24) 

0.20(0.55) 

0.16(0.68) 

0.02(0.09) 

Table  9.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  40  AND 
UC  =  75% 


M.B(M.S.F)  f,r  /?„ 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

-0.06(0.20) 

-0.06(0.  IS) 

-0.1 2(0.  IS) 

-0.05(0.13) 

15  OBS 

-0.007(0.09) 

(— IW 

-0.06(0.06) 

-0.03(0.06) 

30  OBS 

-0.02(0.04) 

-0.01(0.04) 

-0.04(0.03) 

-0.01(0.02) 
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III.  A  WEIBULL  REGRESSION  MODEL 


A.  MODEL 

Suppose  there  are  M  observers.  Each  observer  is  presented  with  N,  targets.  Let  Utj 
be  the  time  it  takes  for  observer  i  to  detect  target  j  .  Let  x,A  ,  xlj7  ,  xl/s ,...,  xIJP  be  the  values 
of  explanatory  variables  which  may  effect  U,,  (  e.g.,  terrain,  atmospheric  conditions,  fa¬ 
tigue  etc.) 

Assume  that  L'v  are  independent  random  variables  with  Weibull  distributions 

P{Uij<  t)  =  \  —  exp[  —  (re-2(/£)e  {'] 

r  -<n  (3-J) 

=  1  -  exp[  -  (tlmj)e  ] 

where  n  =  e'-i  with 

p 

Zij£  =  .  O-2) 

k=\ 

It  follows  from  equation  (3.1)  that  the  distribution  of  the  ln-detection  time  is 
P{  In  L’ij  <  t}  =  P[i)j  <  e'} 

=  1  —  exp{  —  [e'e~  ^  (3.3) 

=  1  -  exp{  -  exp[(r  -  *■]}  . 

The  derivative  of  equation  (  3.3  )  with  respect  to  t  is  following 

P{  In  1'.,  e  dr)  =  exp{  -  exp[(r  -  ■lijB)e~  exp{(r  -  x  S)e~  5,]^_  51 

■"  L  1  _f  (3.4) 

=  expj  -  exp[(r  -  £ijP)e~  J']}  exp{(r  -  £(i§)e~ 51  -  ;,}  . 

When  the  /'*  observer  is  presented  with  his  j*  target,  he  gets  a  length  of  time  called  an 
opportunity  time  0„  to  detect  it.  The  observer  either  successfully  detects  the  target 
within  the  opportunity  time  or  is  unsuccessful.  Data  for  the  observer  consist  of  times 
of  detections  for  the  successes  and  the  lengths  of  opportunity  times  for  the  failures.  For 
each  i  =  1.2,3 . M  ,  j  =  1,2.3 . V,  ,  let 

}]j  =  min(  In  L)..  In  0(;)  (3.5) 
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1  if  Utj  <  O'j 
0  otherwise  . 


(3.6) 


The  YtJ  are  the  censored  ln-detection  times  and  Ay  is  an  indicator  of  whether  or  not  the 
ln-time  to  detect  the  j'h  target  by  the  i*  observer  is  censored.  Let 

S', 

Q  =  (3.7) 

7=1 


be  the  number  of  targets  detected  by  observer  /.  The  next  section  of  this  chapter  will 
discuss  estimation  of  the  parameters  {/?*}  \_Rcf.  4],  These  parameter  estimates  might  be 
used  to  predict  future  times  to  detection  for  an  observer  given  his  past  performance.  In 
the  last  section  of  the  chapter  results  of  a  simulation  study  of  the  estimation  precedure 
will  be  given. 

B.  ESTIMATION 

1.  The  likelihood  Equation  and  Maximum  Likelihood  Estimators. 

The  likelihood  for  the  i'h  observer  is 


UP-  £>  =  1  ~  *J)  exPl  ”  exPL<-l'y  “  -ijP)e  ?']}  •  (-'•&) 

7=1 

The  In  likelihood  for  the  i!h  observer  is 


In  L//?.  c)  = 


£ijP)e  1  -  Q  -  expffvy  -  XijP)e  *'])  . 


(3.9) 


Since  the  observers  are  assumed  to  be  independent,  the  In  likelihood  for  all  observers  is 


V/ 

In  L(§_,  £)  =  Yin  Lt(p,  i) 


i=  1 

.V/  'V, 


=  -  Xfjfle  J'  -  Q  -  exp[(fy  -  £jjP)e  fl]} 


'=!  y=) 


(3.10) 


a.  Newton's  Procedure  to  solve  for  pk 

The  partial  derivative  of  equation  (3.10)  with  respect  to  /?*  is 


M  ' 

jjr-  In  L(/>.  c)  =  Y  ^{A/y[(  -  xijk)e~  "]  -  exp[(j-y-  -  ^)e~  ?'](  -  xiJke~  ’*')] 


i=l  j=\ 


M  -V' 

=  X  ^  i  “  Y +  exP[0'V  ~  • 


(3-11) 


‘=  i  ;=i 


A  Newton  procedure  to  solve  equation  (3.1 1)  for  ^  assuming  {£,}  known  uses  the  second 
derivatives 


-T-l-  =  ^  Y- exp[0v;  -  Xji)e~  *'  •  (3. 1 2j 


<=1  y=l 


vij  =  v  exp[tvv  -  ■£/£)*  *']  ; 


(3.13) 


* ijh  =  xijhe  "wij 


The  Newton  procedure  to  solve  equations  (3.1 1)  for  £  can  be  written  as 
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(3.15) 


0  = 


i=i  j=\  h= 1 
\t  \ 


<=i  ;=i 


which  are  of  the  form  of  the  normal  equations  for  Least  Squares  regression  with  de¬ 
pendent  variables 


:ij  ~ 


and  independent  variables 

uijk  =  wifX‘jke~  • 

b.  Newton  Procedure  to  solve  for  c, 

The  partial  derivative  of  equation  (3.10)  with  respect  to  i,  is 


(3.16) 


(3.17) 


-jr  In  L'ifi,  £,)  =  ^[{A,,[Cvv  -  ’*•(  -  1)  -  1]} 

7=1 

~  {  exp[lliy  “  "Xty  “  &ij §)e~  ?'(  -  1 )}]  (3. 1 8) 

v, 

=  -  Q  +  *'{  -  Ay  +  exp[Cvy  -  £ijfl)e~  !']} 

7=1 

A  Newton  procedure  to  solve  equation  (3.15)  for  l,  asuming  p  known  uses  the  second 
derivative 
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~~T  ln  L‘(P'  ii)  *  -  ~  *Se  S,(  -  AIJ  +  exP[0'V  -  Xyp)e  J']} 

cz-i  7=1 

,v, 

+  -  ±Se~  *'{  exP[C Vij  -  ZijP)e~  ’  JO-'//  -  £ijp)e~  ='(  -  I)} 

7=1 

A' 

=  ~  Q  -  Yoiy  -  ZijfiYe~21  expfCty  -  ^jp)e~  *'] 


(3.19) 


if  jr  In  L,(/J,  £,)  =  O-  Therefore  the  Newton  equation  for  is 


.A . 

0sB7Tln  L/(i-  £)  =  -  ci  +  XA'>'lV  ~  ’"{  -  Ay  +  exp[fvl7  -  £ ijp)e~  f]j 

*'  7=1 

+  (  ~  Ci  ~  -  *7 £)C~  ’f  exP[Cv'i/  -  Xv&)e~  •■]})[{,  -  w] 


(3.20) 


where  c°  is  the  current  value  of  <£..  Solving  results  in  the  equation 


,  a; 

•)  +Ci-  Yj^ij  -  £iji)e~  ?'[  -  Aij  +  exp{Cv(;  -  £ijp)e~  f,]]| 

o  S _  7=1  J 


-  c,  -  V{[0'y  -  Xy£)e  J']2  cxp[(iv  -  £ij§)c  ?’]J 


(3.21) 


2.  Initial  Condition 

In  this  subsection,  we  describe  a  rough  way  to  provide  initial  estimates  ofx/? 
to  start  the  iterative  Newton  procedure  of  the  previous  section.  It  follows  from  equation 
(3.1)  that 
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=  l-exp[-(/eV^)f'<] 

=  1  —  exp[  — 


(3.22) 


Hence, 


£[(!••/’ ‘I- 


(3.23) 


and 


In  £[(C/"{']=^-f' 

As  a  result,  for  all  the  observed  i]  censored  or  uncensored,  we  will  put 


(3.24) 


We  will  take  <f? 
u  is 


0  for  convenience.  Thus  the  initial  value  of  x./?°  for  each  observation 


&ijP'  =  In  u( 

3.  Recursive  Procedure  to  find  {/?*)  and 

1.  Put  T  =  0  for  /  =  1.2.3 . M 

2.  Compute  the  initial  value 

.Sjjf  =  e"  In  Vjj 

3.  Iteration 

a.  Compute  dependent  and  independent  variables  for  regression 

ivij  =  \  exp{Civ  -  £<£)€~  5,1 

uijk  =  **  )y  xijk  e~  "  ( dependent  variable) 

-A  a+Wt: 


(3.261 


(3.27) 


(3.28) 


— - t-  IVjjjCj-P  {independent  variable) 


Let 
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p  1 

Ml 

“in 

j: 

to 

ull/> 

w121 

w122 

w12p 

Zl-'i 

221 

U1.V,2 

W1  A> 

u211 

w212 

w21p 

Zi  v 

L  = 

♦ 

* 

4.  .  »*> 

“2\2I 

w2  ,Vj2 

W2.\> 

2  Ml 

UM11 

WM12 

uM1/> 

;.irv„ 

“V.VM1 

MW\w2 

U M  Vwf 

b.  Compute  regression  estimates 

p  =  (ltd~'ltz 


c.  Recompute  the  initial  value 


Xijfi  = 


d.  Cpdate  :  as  m  equation  (3.21) 


r  v  1 

*  +  C,  -  Yiy,  -  Xi£)c~  *'•[  -  \j  +  exp {(.!•//  -  ‘  (]  • 


t  _  S'-  _ 
S-  «s-  — 


l 


-  Q  -  V { [(ty  -  XijP)c  *'}2  exp[0'iy  -  £ijH)e  *'] f 


y=i 


e.  Put 


=  f  / 


f.  Return  to  step  3 


4.  Iterate  until  the  following  conditions  are  satisfied 


f  I  ~  ft  I 
\ 


<  l£  —  4.  k  =  l,2.3,...,/> 


(3.29) 


(3.30) 


T* 


max 


C.  SIMULATION  PROCEDURE  AND  RESULTS 

1.  Simulation 

The  number  of  replications  for  each  simulation  experiment  is  100.  Each  simulation  ex¬ 
periment  of  100  replications  starts  with  the  same  random  number  seed.  The  mean  bias.its  standard 
error  and  mean  square  error,  standard  error  of  mean  square  error  from  R  =  100  replications  are 
computed  as  in  equations  (2.49a)  through  (2.50b).  In  each  of  the  experiments  the  number  of 
covariates  for  each  observer  is  2. 

The  simulation  experiment  to  study  the  sampling  properties  of  the  estimators  of  {/?*}  and 
is  as  follows. 

1.  Give  arbitrarv  constant  values  for  the  true  values  of  /?A  and  £..  In  all  the  simulations  /?,  = 
0.2.  p2  =  0.3  and  =  0. 

2.  The  same  opportunity  time,  0„.  i  =  1.2,3 . M.  and  j  =  1,2, 3,. ...A',,  is  used  for  \1  observers. 

The  values  of  the  constant  opportunity  time  are  1.3.  2.5,  4.0. 

3.  Generate  random  numbers 

a.  Generate  independent  normal  random  numbers  with  mean  1  and  variance  0.5.  x  •  .  i  = 

1.2.3 M ../  =  1.2.3 . A',. 

b.  Generate  independent  normal  random  numbers  with  mean  2  and  variance  1  ,  /  = 

1.2.3 M.y  =  1.2,3 . V. 

c.  Generate  exponential  random  numbers  with  mean  1,  IT,,,  i  =  1.2.3 . Nl.  j  =  1.2.3 . 

A'.. 

4.  Compute  ,u  .  as  follows  : 

=  exp(j3 j X:j i  +  p2xijZ) .  (3.3 1  i 

5.  Compute  the  detection  time  that  it  takes  for  observer  i  to  detect  target  j  as 

=  .  (3.32) 

6.  Compare  In  l',,  with  In  O  .  then  choose  smaller  one  for  the  data 

Yji  =  min(  In  L In  Or  ) 

and  compute  A  as  equation  (3.7) 

7.  Compute  the  initial  value  as  in  equation  (3.27) 

S.  Compute  the  values  of  the  regression  variables  .  zr  and  uri  as  in  equation  (3.16)  and  (3.17) 

9.  Compute  regression  estimates  as  in  equation  (3.29) 

10.  Recompute  the  initial  value  with  new  /?  value  as  in  equation  (3.30) 

1 1.  Compute  g  as  in  equation  (3.21) 

12.  Iterate  the  procedure  until  the  differences  between  successive  values  of  /?, ,  (32  and  T  are  small. 
(  less  than  10  4  ). 

2.  Results 

In  this  section  results  from  the  simulation  experiments  will  be  reported.  For  R  repli¬ 
cations.  the  uncensonng  level  (  UC  )  for  an  experiment  is  computed  as  follows  : 
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LC  = 


R  M 


r=l  i=l 


/?  U 


V  V 


r=l  /=! 


where  C,{r )  is  the  number  of  the  targets  detected  by  obsen'er  /  in  replication  r  and  S,{r)  is  the 
number  of  targets  presented  to  observer  i  in  replication  r.  The  UC's  for  the  simulations  using  the 
same  opportunity  time  are  then  averaged  to  obtain  the  mean  UC  for  that  opportunity  time.  Figures 
19  through  45  show  histograms  of  (  /?*  -/?*),  k  =  1,2  and  (  -  £  ).  Tables  10  through  21  show- 

mean  square  errors  and  mean  biases  for  each  of  the  estimates,  their  standard  errors  appear  in  pa¬ 
renthesis  below.  Tables  10.  14  and  18  present  all  the  simulation  results  for  different  values  of  op¬ 
portunity  times.  The  other  tables  present  the  same  results  in  a  more  convenient  fashion.  The  mean 
biases  are  displayed  with  the  mean  square  error  in  parenthesis;  for  each  table,  the  opportunuty 
time(O)  and  average  uncensoring(UC)  for  that  are  given  at  the  top  of  table. 

•  Results  for  /?,  and  p2 

The  histograms  for  /?,  -  /?,  and  (]2  -  jS2  tend  to  be  somewhat  centered  around  0.  For 
small  numbers  of  targets  and  observers  the  histograms  tend  to  be  slightly  skewed  to  the  left. 
Increasing  the  observation  time  has  little  effect  on  the  mean  bias  and  mean  square  error.  In¬ 
creasing  the  number  of  observers  for  a  fixed  number  of  targets  tends  to  decrease  the  mean 
square  error  bui  has  less  effect  on  the  mean  bias.  Increasing  the  number  of  targets  for  a  fixed 
number  of  observers  decreases  the  mean  bias  and  mean  square  error.  The  standard  errors  of 
the  mean  biases  mean  square  errors  are  large. 

•  Results  for 

The  histograms  of  c  -  c  tend  to  be  centered  about  0.  There  is  a  tendency  for  slight 
skewness  to  the  right  for  small  numbers  of  observers  and  targets.  Changing  the  opportunity 
time  has  little  effect  on  the  mean  bias  and  mean  square  error.  The  more  targets  there  are  for 
observer  i.  the  smaller  the  mean  bias  and  mean  square  error  for  The  standard  errors  of  the 
mean  biases  and  mean  square  errors  are  large. 
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Table  10.  MEAN  BIAS,  MEAN  SQUARE  ERROR  AND  STANDARD  ERROR  AT 
0=  1.3  AND  UC  =  45% 


Number  of 
Observers 

Number  of 
Targets 

@2 

MB 

(S.E) 

M.S.E 

(S.E) 

MB 

(S.E) 

M.S.E 

(S.E) 

MB 

(S.E) 

M.S.E 

(S.E) 

5 

5 

-0.06 

(0.03) 

0.07 

(0.01) 

-0.02 

(0.02) 

0.03 

(0.01) 

-0.05 

(0.01) 

0.02 

(0.01) 

15 

-0.02 

(0.02) 

0.02 

(0.01) 

-0.02 

(0.01) 

0.01 

(0.001) 

-0.02 

(0.01) 

0.01 

(0.001) 

30 

-0.003 

(0.01) 

0.02 

(0.01) 

-0.02 

(0.01) 

0.01 

(0.002) 

-0.03 

(0.01) 

0.01 

(0.001) 

50 

-0.01 

(0.001) 

0.01 

(0.001) 

-0.01 

(0.01) 

0.003 

(0.001) 

-0.01 

(0.01) 

15 

5 

-0.05 

(0.01) 

0.02 
(0  003) 

-0.04 

(0.01) 

0.01 

(0.001) 

-0.02 

(0.01) 

0.01 

(0.001) 

15 

-0.03 

(0.01) 

0.01 

(0.001) 

-0.02 

(0.01) 

0.003 

(0.001) 

-0.01 

(0.01) 

0.002 

(0.0003) 

30 

-0.01 

(0.01) 

0.01 

(0.001) 

-0.01 

(0.004) 

0.002 

(0.001) 

-0.01 

(0.004) 

0.002 

(0.0002) 

50 

-0.01 

(0.001) 

0.003 

(0.0004) 

-0.01 

(0.003) 

0.001 

(0.0004) 

-0.01 

(0.004) 

0.001 

(0.0002) 

30 

5 

-0.05 

(0.01) 

0.01 

(0.002) 

-0.04 

(0.01) 

0.005 

(0.001) 

-0.02 

(0.005) 

0.003 

(0.001) 

15 

-0.02 

(0.01) 

0.003 

(0.0004) 

■  ■ 

0.001 

(0.0002) 

-0.01 

(0.003) 

0.001 

(0.0002) 

30 

-0.01 

(0.01) 

0.003 
(0  1)004) 

-0.01 

(0.003) 

0.001 

(0.0001) 

-0.004 

(0.003) 

0.001 
(O.OOd  1 1 

50 

-0.01 

(0.004) 

0.002 

(0.0002) 

-0.01 

(0.003) 

0.001 

(0.0001) 

-0.01 

(0.002) 

0.001 
(O.OdOl ) 
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Table  11.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0=1.3  AND 
UC  =  45% 


.\f.B(\l.S.E)  for  /?, 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

50BS 

-0.04(0.06) 

-0.02(0.02) 

-0.004(0.02) 

-0.003(0.01) 

150BS 

-0.05(0.02) 

-0.02(0.007) 

-0.01(0.005) 

-0.01(0.003) 

300  BS 

-0.008(0.01) 

-0.05(0.003) 

-0.006(0.003) 

-0.01(0.002) 

Table  12.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0=1.3  AND 
UC  =  45% 


M.B(\f.S.E)  for 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

-0.03(0.02) 

-0.02(0.00") 

-0.02(0.005) 

-0.01(0.002) 

15  OBS 

.0.04(0.007) 

-0.009(0.003 

-0.01(0.002) 

-0.01(0.001) 

30  OBS 

.0.(14(0.004) 

-0.01(0.001) 

-0.01(0.001) 

-0.009(0.0006) 

Table  13.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0=1.3  AND 
UC  =  45% 


M.BiM.S.E)  for  l 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

-0.05(0.14) 

-0.02(0.007) 

-0.02(0.005) 

-0.006(0.004) 

15  OBS 

-0.01(0.005) 

-0.009(0.002) 

-0.005(0.002) 

-0.008(0.001) 

30  OBS 

-0.02(0.002) 

-0.01(0.001) 

-0.004(0.0007) 

0.008(0.0006) 
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Table  14.  MEAN  BIAS,  MEAN  SQUARE  ERROR  AND  STANDARD  ERROR  AT 
O  =  2.5  AND  UC  =  75% 


Number  of 
Observers 

Number  of 
Targets 

k 

k 

i 

MB 

(S.E) 

M.S.E 

(S.E) 

M.B 

(S.E) 

M.S.E 

(S.E) 

M.B 

(S.E) 

M.S.E 

(S.E) 

5 

5 

-0.04 

(0.02) 

0.06 

(0.01) 

-0.03 

(0.01) 

0.02 

(0.003) 

-0.05 

(0.01) 

0.14 

(0.01) 

15 

-0.02 

(0.02) 

0.02 

(0.003) 

-0.02 

(0.01) 

0.01 

(0.001) 

-0.02 

(0.01) 

0.01 

(0.001) 

30 

-0.004 

(0.01) 

0.02 

(0.003) 

-0.02 

(0.01) 

mw 

-0.02 

(0.01) 

0.01 

(0.001) 

50 

-0.003 

(0.01) 

0.01 

(0.001) 

-0.01 

(0.01) 

0.002 

(0.0003) 

-0.01 

(0.01) 

0.004 

(0.001) 

15 

5 

-0.05 

(0.01) 

0.02 

(0.003) 

-0.04 

(0.01) 

0.01 

(0.001) 

-0.01 

(0.01) 

0.005 

(0.001) 

15 

-0.02 

(0.01) 

0.01 

(0.001) 

-0.02 

(0.01) 

0.003 

(0.0004) 

-0.01 

(0.01) 

0.002 

(0.0003) 

30 

-0.01 

(0.01) 

0.01 

(0.001) 

-0.01 

(0.004) 

0.002 

(0.0002) 

-0.01 

(0.004) 

0.002 

(0.0002) 

50 

-0.01 

(0.01) 

0.003 

(0.0004) 

-0.01 

(0.003) 

0.001 

(0.0002) 

-0.01 

(0.003) 

0.001 

(0.0002) 

30 

5 

-0.01 

(0.01) 

0.01 

(0.001) 

-0.04 

(0.01) 

0.004 

(0.001) 

-0.02 

(0.01) 

0.003 

(0.0004) 

15 

-0.05 

(0.01) 

0.003 

(0.0004) 

-0.02 

(0.003) 

0.001 

(0.0002) 

-0.01 

(0.003) 

0.001 

(0.0002) 

30 

-0.01 

(0.01) 

0.003 

(0.0003) 

-0.01 

(0.003) 

0.001 

(0.0001) 

-0.004 

(0.003) 

0.001 

(0.0001) 

50 

-0.01 

(0.004) 

0.0012 

(0.0002) 

-0.01 

(0.002) 

0.001 

(0.0001) 

-0.01 

(0.002) 

0.001 

(0.0001) 
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Table  15.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  2.5  AND 
UC  =  75% 


M.B(M.S.E)  for  /?, 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

-0.06(0.07) 

-0.02(0.20) 

-0.003(0.02) 

-0.005(0.01) 

15  OBS 

-0.05(0.02) 

-0.03(0.007) 

-0.01(0.006) 

-0.008(0.003) 

30  OBS 

-0.05(0.01) 

-0.02(0.003) 

-0.006(0.003) 

-0.01(0.002) 

Table  16.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  2.5  AND 

UC  =  75% 


\1.B\  M.S.E)  for  j}2 

5  TGT 

15  TGT 

30  TGT 

50  IGF 

5  OBS 

-0.02(0.03) 

-0.02(0.007) 

-0.02(0.005) 

-0.006(0.003) 

15  OBS 

-0.04(0.007) 

-0.02(0.003) 

-0.0 1(0.002) 

-0.01(0.001) 

30  OBS 

-0.04(0.004) 

-0.02(0.001) 

-0.01(0.001) 

-0.006(0.001) 

Table  17.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  2.5  AND 
UC  =  75% 


M.B(M.S.E)  for  l 

5  TGT 

15  TGT 

30  TGT 

50  TGT 

5  OBS 

-0.046(0.02) 

-0.03(0.008) 

-0.03(0.004) 

-0.006(0.004) 

15  OBS 

-0.02(0.01) 

-0.01(0.002) 

-0.006(0.002) 

-0.009(0.001) 

30  OBS 

-0.02(0.003) 

-0.01(0.001) 

-0.004(0.0008) 

-0.009(0.0006) 
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Table  18.  MEAN  BIAS,  MEAN  SQUARE  ERROR  AND  STANDARD  ERROR  AT 
0  =  4.0  AND  UC  =  89% 


Number  of 

Number  of 

h 

l 

Observers 

Targets 

MB 

(S.E) 

Nl.S.E 

(S.E) 

MB 

(S.E) 

M.S.E 

(S.E) 

M.B 

(S.E) 

M.S.E 

(S.E) 

5 

-0.03 

(0.01) 

0.02 

(0.003) 

-0.03 

(0.02) 

0.03 

(0.01) 

-0.07 

(0.01) 

0.02 

(0.003) 

< 

15 

-0.01 

(0.01) 

0.01 

(0.001) 

-0.02 

(0.01) 

0.01 

(0.001) 

-0.02 

(0.01) 

0.01 

(0.001) 

J 

30 

-0.01 

(0.01) 

0.01 

(0.001) 

-0.02 

(0.01) 

0.01 

(0.001) 

-0.03 

(0.01) 

11 

50 

-0.002 

(0.01) 

0.01 

(0.001) 

-0.01 

(0.01) 

0.004 

(0.001) 

-0.01 

(0.01) 

0.003 

(0.0004) 

5 

-0.03 

(0.01) 

0.01 

(0.001) 

-0.04 

(0.01) 

0.01 

(0.001) 

-0.04 

(0.01) 

0.01 

(0.001) 

15 

15 

-0.02 

(0.01) 

0.003 

(0.0004) 

-0.02 

(0.001) 

0.004 

(0.0001) 

-0.02 

(0.0004) 

0.002 

(0.0003) 

30 

-0.01 

(0.001) 

0.002 

(0.0003) 

-0.01 

(0.001) 

0.002 

(0.0003) 

-0.01 

(0.004) 

0.001 

(0.0002) 

50 

-0.01 

(0.004) 

0.002 

(0.0002) 

-0.01 

(0.04) 

0.002 

(0.0002) 

-0.01 

(0.003) 

0.001 

(0.0002) 

5 

-0.03 

(0.01) 

0.004 

(0.001) 

-0.04 

(0.01) 

0.006 

(0.001) 

-0.04 

(0.01) 

0.004 

(0.001) 

30 

15 

-0.01 

(0.004) 

0.002 

(0.0002) 

-0.02 

(0.004) 

0.002 

(0.0003) 

-0.02 

(0.003) 

0.001 

(0.0002) 

30 

-0.03 

(0.0O03) 

0.001 

(0.00) 

-0.01 

(0.0004) 

0.001 

(0.00) 

-0.01 

(0.0003) 

0.001 

(0.00) 

50 

-0.01 

(0.0003) 

0.001 

(0.0U) 

-0.01 

(0.003) 

0.001 

(0.0001) 

-0.01 

(0.002) 

0.001 
(0.000 1 ) 

Table  19.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  4.0  AND 
UC=  89% 


Table  20.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  4.0  AND 
VC  =  89% 


Table  21.  MEAN  BIAS  AND  MEAN  SQUARE  ERROR  AT  0  =  4.0  AND 
UC  =  89% 
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IV.  CONCLUSION 


This  thesis  considers  two  models  for  the  times  until  detection  of  targets.  Each  model  has  M 
observ  ers.  The  i'h  observer  is  presented  with  A',  targets.  In  the  model  of  Chapter  2,  observer  i  has  a 
random  variable  A,  which  reflects  his  ability  to  detect  a  target.  The  random  variables  {A,}  are  inde¬ 
pendent  identically  distributed  having  a  gamma  distribution.  Given  A„  the  times  to  target  detections 
for  observer  i  are  conditionally  independent  Weibull  random  variables  with  known  parameters. 
Simulation  experiments  indicate  that  increasing  the  number  of  observers  for  a  fixed  number  of  tar¬ 
gets  provides  the  greatest  decrease  in  the  mean  bias  and  mean  square  error  of  the  estimates  of  the 
parameters  of  the  gamma  distribution  that  describes  the  variation  between  individuals.  This  is  not 
surprising,  since  observing  more  individuals  sampled  from  a  fixed  population  should  better  estimate 
properties  of  that  population.  The  model  of  Chapter  3  is  a  Weibull  regression  model.  In  this  case 
the  simulation  experiments  indicate  that  increasing  the  number  of  targets  for  a  fixed  number  of 
observers  provides  the  greatest  decrease  in  the  mean  bias  and  mean  square  error  of  the  estimators. 

A  topic  for  future  investigation  is  to  combine  the  two  models  and  estimation  procedures  to 
provide  estimates  for  a  hierarchical  gamma  Weibull  regression  model.  Another  topic  is  to  investi¬ 
gate  using  the  fitted  hierarchical  model  to  predict  future  performance  of  the  observers. 
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FREQUENCY  FREQUENCY 


APPENDIX  A.  HISTOGRAMS  FOR  THE  ESTIMATED  GAMMA 

PARAMETERS 


5  OBSERVERS  (  UC  :  37  PERCENT  ) 

5  TARGETS  15  TARGETS 


30  TARGETS 


ETTA  -  A  TRUE 
50  TARGETS 


ETA  -  ETA  TRUE 


Figure  1.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=10, 

UC  =  37%):  5,15,30,50  targets  for  5  observers  (fj  -  t]) 
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FREQUENCY  FREQUENCY 


5  OBSERVERS  (  UC  : 

5  TARGETS 


62  PERCENT  ) 

15  TARGETS 


ETA  -  ETA  TRUE 
30  TARGETS 


ETA  -  ETA  TRUE 


ETA  -  ETA  TRUE 
50  TARGETS 


ETA  -  ETA  TRUE 


Figure  2.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  25, 

UC  =  62%):  5,15,30,50  targets  for  5  observers  (tj  —  rj) 
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FREQUENCY  FREQUENCY 


Figure  3.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  40, 

UC  =  75%):  5,15,30,50  targets  for  5  observers  {rj  —  tj) 
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FREQUENCY  FREQUENCY 


15  OBSERVERS  (  UC  :  37  PERCENT  ) 


5  TARGETS 


ETA  -  ETA  TRUE 
30  TARGETS 


ETA  -  ETA  TRUE 


15  TARGFT3 


ETA  -  ETA  TRUE 
50  TARGETS 


DA  -  ETA  TRUE 


Figure  4.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=10, 

UC  =  37%):  5,15,30,50  targets  for  15  observers  (fj  —  tj) 
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frequency  frequency 


1 


5  TARGETS 


/ 

15  TARGETS 


30  TARGETS 


HA  -  ETA  TRUE 
50  TARGETS 


ETA  -  ETA  TRUE  ETA  -  HA  TRUE 


Figure  5.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  25, 

UC  =  62%):  5,15,30,50  targets  for  15  observers  (tj  —  ?/) 
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FREQUENCY  FREQUENCY 


Figure  6.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  40, 

UC  =  75%):  5,15,30,50  targets  for  15  observers  (>/  —  r\) 


m 
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FREQUENCY 


30  OBSERVERS  (  UC  :  37  PERCENT  ) 

5  TARGETS  15  TARGETS 


Figure  7.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=10, 

UC  =  37%):  5,15,30,50  targets  for  30  observers  (rj  -  tj) 


FREQUENCY  FREQUENCY 


30 


Figure  8.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  25, 

UC  =  62%):  5,15,30,50  targets  for  30  observers  {rj  -  tj) 
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FREQUENCY  FREQUENCY 


30  OBSERVERS  (  UC 

5  TARGETS 


75  PERCENT  ) 

15  TARGETS 


30  TARGETS 


ETA  -  ETA  TRUE 
50  TARGETS 


ETA  -  ETA  TRUE  ETA  -  ETA  TRUE 


Figure  9.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  40, 

UC  =  75%):  5, 1 5,30,50  targets  for  30  observers  -  tj) 
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frequency  frequency 


5  OBSERVERS  (  UC  :  37  PERCENT  ) 


5  TARGETS 


15  TARGETS 


BETA  -  BETA  TRUE 


BETA  -  BETA  TRUE 


Figure  10.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=10, 

UC=  37%):  5,15,30,50  targets  for  5  observers  (/?„  —  /?„) 
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FREQUENCY  FREQUENCY 


5  OBSERVERS  (  UC 

5  TARGETS 


62  PERCENT  ) 

15  TARGETS 


BETA  -  BETA  TRUE 
50  TARGETS 


50  TARGETS 


BETA  -  BETA  TRUE  BETA  -  BETA  TRUE 


Figure  11.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  25, 

UC  =  62%);  5,15,30,50  targets  for  5  observers  (j?0  -  /?„) 
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FREQUENCY  FREQUENCY 


5  OBSERVERS  ( 

5  TARGETS 


BETA  -  BETA  TRUE 
30  TARGETS 


UC  :  75  PERCENT  ) 

15  TARGETS 


-2  0  2 
BETA  -  BETA  TRUE 

50  TARGETS 


BETA  -  BETA  TRUE  BOA  -  BETA  TRUE 


Figure  12.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  40, 

UC  =  75%):  5,15,30,50  targets  for  5  observers  (/?„  -  /?„) 


50 


Figure  13.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=10, 

UC  =  37%):  5,15,30,50  targets  for  15  observers  (/?„  —  /?„) 
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FREQUENCY  FREQUENCY 
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Z3 

S 

£ 


XL 


s 
£  . 


j_d 


tL 


-2 


-2 


BETA  -  BETA  TRUE 


BETA  -  BETA  TRUE 


Figure  15.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  40, 

UC  =  7 5 % ) :  5,15,30.50  targets  for  15  observers  (/?0-/?0) 


53 


IMS 


5  TARGETS 


15  TARGETS 


BETA  -  BETA  TRUE 
30  TARGETS 


BETA  -  BETA  TRUE 


50  TARGETS 


T 

J 


BETA  -  BETA  TRUE 


b 

z 

UJ 

z> 

o 

u 

e 


BETA  -  BETA  TRUE 


Figure  16.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=10, 

UC  =  37%);  5,15,30,50  targets  for  30  observers  (/?„  -  /?„) 


-I _ I— mUj 

-2  0 
BETA 


h  ■ 

BETA  TRUE 


BETA  -  BETA  TRUE 


Figure  17.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  25, 

UC  =  62%):  5, 1 5,30,50  targets  for  30  observers  (/?0  —  /?0) 


FREQUENCY  FREQUENCY 


APPENDIX  B.  HISTOGRAM  FOR  THE  WEIBULL  REGRESSION 

PARAMETERS 


5  OBSERVERS  (  UC  :  45  PERCENT  ) 


5  TARGETT5 


30  TARGETS 


15  TARGETS 


50  TARGETS 


Figure  19.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=1.3. 

UC  =  45%):  5,  15,  30,  50  targets  For  5  observers  (/?,  -  /?,) 


FREQUENCY  FREQUENCY 


5  OBSERVERS  (  UC 

5  TARGETS 


:  75  PERCENT  ) 

15  TARGETS 


BETA1  -  BETA!  TRUE 


30TARGETS 


8ETA1  -  BETA!  TRUE 


BETA1  -  BETA1  TRUE 
50TARGETS 


Figure  20.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  2.5, 

UC  =  75%);  5,  15,  30,  50  targets  for  5  observers  (/?,  —  /?,) 


58 


5  OBSERVERS  ( 


5  TARGETS 


BETA)  -  BETA)  TRUE 
30  TARGETS 


UC  :  89  PERCENT  ) 


15  TARGETS 


50  TARGETS 


-0.4 


-0.2  0  0.2 
BETA1  -  8ETA1  TRUE 


0.4 


Figure  21.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0-4.0, 

UC  =  89%):  5,  15,  30,  50  targets  for  5  observers  (/?,  ~  Pi) 
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FREQUENCY  frequency 
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FREQUENCY  frequency 


61 


FREQUENCY  FREQUENCY 


FREQUENCY  FREQUENCY 


n_i_ 


-0.8 


-0.4  0 

BE7A1  -  BETA!  TRUE 


0.4 


_ I 

0.8 


ll _ L_ 


-0.6 


-0.4  0 

BETA1  -  BETA1  TRUE 


Figure  26.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  2.3, 

UC  =  75%):  5,  15,  30,  50  targets  for  30  observers  (/?,  -  /?,) 
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i _ i _ i _ jimmiii-,  ,, 

-0.4  -0.2  0  0.2  04 

BETA1  -  BETA1  TRUE 


3 

a 

UJ 


e 


o 


_ I _ L 

-0.4 


-0 


2  0 
BETA1  -  BC7AI  TRUE 


0.2 


1 


Figure  27.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  4.0, 

UC  =  89%):  5,  15,  30,  50  targets  for  30  observers  (j?,  — /?,) 


0  4 
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i 


u\ 


5  TARGETS 


15  TARGETS 


-0.4 


0.4 


BETA2  -  BETA2  TRUE 
30  TARGETS 


-I _ L 


-0.4 


tk. 


0.4 


BETA2  -  BETA2  TRUE 


z; 

LJ 

o 


-0.4 


BETA2  -  BETA2  TRUE 


0.4 


Figure  28.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=1.3, 

UC  =  45%):  5,  15,  30,  50  targets  for  5  observers  {p2-  /?2) 
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Figure  29.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  2.5, 

UC  =  75%):  5.  15,  30,  50  targets  Tor  5  observers  (jj2  -  fl2 ) 


I 


frequency  frequency 


5  OBSERVERS  (  UC 


89  PERCENT  ) 


5  TARGETS 


15  TARGETS 


BETA2  -  BETA2  TRUE 


30  TARGETS 


50  TARGETS 


BETA2  -  BETA2  TRUE 


Figure  30.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  4.0, 

UC  =  89 %):  5,  15,  30,  50  targets  for  5  observers  (Pi  ~  Pi) 
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FRCQUENCY  rREQUENCT 


0  0  2 
enA2  -  BtTA2  true: 


0.4 


-0  4 


“0  2  0  0  2 
BFTA2  -  BETA 2  TRUE 


COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  4.0, 

CC  =  89%):  5,  15,  30,  50  targets  for  15  observers  if>2  —  fi2 ) 


frequency  frequency 


o.a 


-0.4 


-0.8 


0  4 


0  8 


BETTA2  -  BETA2  TRUE 


-0.4  0  0  4  Ct' 

BETA2  -  8ETA2  TRUE 


Figure  35.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  2.5, 

UC  =  75%):  5,  15,  30,  50  targets  for  30  observers  (P2  -  /?,) 


j - 1 - 1 _ mill  I  IJLIi _ i _ i  ...i,  i  •  ,i 

04  -0.2  0  0.2  0.4  0.8 

BETA2  -  BETA2  TRUE 


—,i  .  .i _ l _ irfllllllh  i  i  i _ i _ i _ i 

-0.4  -02  0  0.2  0.4  0  0 

BE7A2  -  BE1A2  TRUE 


Figure  36.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  4.0, 

UC  =  89%):  5,  15,  30,  50  targets  for  30  observers  (/?2  -  fl:) 
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FRECUENCY 


Figure  37.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=1.3, 

UC  =  45%):  5,  15,  30,  50  targets  for  5  observers  (£,  -  £,) 
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FREQUENCY  FREQUENCY 


5  OBSERVERS  (  UC  :  75  PERCENT  ) 

5  TARGETS  15  TARGETS 


Figure  38.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  2.5. 

UC  =  75%):  5,  15,  30,  50  targets  for  5  observers  (<*,  —  <?,) 
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FREQUENCY  FREQUENCY 


5  OBSERVERS  (  IJC  :  89  PERCENT  ) 

5  TARGETS  15  TARGETS 


-0.4  -0.2  0  0.2  0.4  -0.4  -0.2  0  0  2  0  4 

THAI  -  THAI  TRUE  THAI  -  THAI  TRUE 


Figure  39.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  4.0, 

UC  =  89%):  5,  15,  3U,  50  targets  Tor  5  observers  (£,  -  £,) 
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FREQUENCY  FREQUENCY 


15  OBSERVERS  (  UC  :  45  PERCENT  ) 


5  TARGETS 


15  TARGETS 


50  TARGETS 


Figure  40.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0=1.3, 

UC  =  45%):  5,  15,  30,  50  targets  for  15  observers  (£,—  £,) 
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FRECUENCT  FREQUENCY 


15  OBSERVERS  (  UC  :  75  PERCENT  ) 


5  TARGETS 


1 5  TARGETS 


50  TARGETS 


Sr 


THAI  -  THAI  TRUE 


Figure  41.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  2.5, 

UC  =  75%):  5,  15.  30,  50  targets  for  15  observers  (£,  —  £,) 
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FREQUENCY 


THAI  -  THAI  TRUE  THAI  -  THAI  TRUE 

Figure  42.  COMPARISON  BETWEEN  DIFFERENT  ESTIMATES  (0  =  4.0, 

UC  =  89%):  5,  15,  30,  50  targets  for  15  observers  (£,  —  £,) 


SO 


FREQUENCY  FREQUENCY 
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FREQUENCY  FREQUENCY 
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APPENDIX  C.  SIMULATION  PROGRAM  FOR  THE  ESTIMATION 

OF  GAMMA  PARAMETERS 


VflAE[D]V 

V  R+M  BAE  N ; MU ;TAiH;E;L; OBS iW;Gl;G;YiC; DEL ; 5 ; 
MllM2;NINiL;V;FiAiBiINiJN-,Cl 

[1]  MU+ 4.2 

[2]  TA+ 0 

C3]  IN+NIN+JN+ 2p0 

[4]  0BS<-4O 

[5]  W<r(M,N)p(M*N)  EXPRAND  1 

[6]  Gl+*s>(N  ,M)pG+M  GAMRAND  1.5  0.154 

[7]  Y+<S>MU*  (W*G1  )*  (*TA  ) 

[8]  Cl<r+/C<r+/DEL+YZ®OBS 
C9]  Y<-Yl ®OBS 

[10]  S<r+/Sl<-*H<-(Y-®MU)x*~lxTA 

[11]  Ml+(.  +  /  +  /DEL*H')*(L++/C) 

[12]  M2+(+/+/DEL*(H*2))iL 

[13]  V«-(M2-(Ml*2))-((ol)*2)*6 

[14]  IN  [  2  ]  •«-  (  1*M1  )-0 . 5772 

[15]  ->A2xi(7>o) 

[16]  V+l 

[17]  A2  :  INll]  +  Clx®V')~INL2l 

[19]  LsClN+NIN 

[20]  L7  :NINl2]+INl2}4-(~lx  +  /  (A-C)*B)*(+/ (A<r(*INl2']  )xS)* 
(B+S+*IN(  1]  )  ) 

[21]  JAN-IW[1]  ,JVJW[2] 

[22]  F+C  SUM  JN 


[23]  iVJtf[l]«-Itf[i]  +  (“ix  +  /(.nV[l]-(®B))+F[il]  )* 
(  +  /(l-*IiV[l]*B)-F[;2]  ) 

[24]  ( ( | (iW[l] -NINlll )*INL 1] )>0 . 0001 )a 

(  (  I  (.INI  2]  [2]  )*JiV[2]  )>0 . 0001  ) 

[2  5]  OWIW 

[26]  R+(NIN+(~  1.872  1 . 466 5  )  )  .  (Cl*Mx/V ) 


V 

VSC/M[D]V 

V  B<-C  SOT  INiDD;BB;AA;Dl;D2;I 

[1]  ZJ0«-iO 

[2]  I<-1 

[3]  L3:+£4xt (C[J]=0) 

[4]  Cl*  +  /*U£«-*i7V[l]+I/V[2]  )+BB«-(iC[7]  )-l 

[5]  D2++/AA* (AA+BB)*2 

[6]  +L 5 

[7]  £4:£l«-£>2«-0 

[8]  £5 :DD+DD,D1 ,D2 

[9]  +L3x\(pC)2I+I+l 

[10]  R+DD*-(  (pC)  ,2  )pDD 


V 

V  THESIS  zm 

V  R+K  THESIS  J;I 

[1]  0HL<e466  801747 

[2]  Mt-l+J 

[3]  N<r~l+J 

[4]  1+ l 

[5] 

£6]  RES+UC,  3)p0 

[7]  LO'.RESLIO+M  BAE 

[8]  +L0*\K2I+I+1 
C  93  R+RES 

V 


iRESiMiN 


N 
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APPENDIX  D.  SIMULATION  PROGRAM  FOR  THE  ESTIMATION 
OF  REGRESSION  PARAMETERS 


[1] 

C2] 

[3] 

[4] 

[5] 

C6] 

[7] 

[8] 

[9] 

CIO] 

[11] 

[12] 

[13] 

[14]  LL :  B0<-BN 

[15]  TO+T 

[16]  T2+T1 

[17]  Wl+(*(<y-iy)x*(-ix2'2)))*0.5 

[18]  Ml+WlxXlx*ClxF2) 

[19]  M2+WlxX2x*ClxT2) 

[20]  MM(2,(MxW))p(,Ml),(,M2) 

[21]  Z+Sl,Z<-(WlxIV')+(.ClxDEL)+Wl*2')*Wl 
[2  2]  BN+ ( 1  ( §MU )  + . xMU ) + , x ( ( §MU )  + .  x  Z  ) 

[23]  IV+(XlxBNlll)+(X2xBNL21) 

[24]  RR+(Y-IV)x*Cl*T2) 

[2  5]  Tl+$(N  ,M)pTl<rT<-,T+T+(C-  +  /RRx  (  ClxDEL  )+*RR)  )* 
( (~lxC)-+/ (RR*2 ) x *RR ) 

[26]  A<-(  T  /  |  (BN-BO  )*BN)>0 . 0001 

[27]  S*| (T0-T)*T 

[28]  +LEx\2QQ<,j+j+\ 

[29]  +LLx\ Qv(r/S)>0.0001  ) 

[30]  LE:R+(BN+BN-  0.2  0 . 3 ) , (Cl* (.MxN) ) , ( (+/T)*M) 

V 

V THESIS [0] 

y  THESIS  J;I;RES-,M;N 

[1]  -<-466  801747 

[2]  M+l+J 

[3] 

[4]  1+ 1 

[5]  RES<-(K,  4)p0 

[6]  L0  : RES II ;  ]-«-W  S71E2  /V 

[7]  -*-L0x  \K^I-*-I+l 

[8] 


VBAE2 [□] y 

y  R+M  BAE 2  N;OBSiBOiBNiWiUliMl;M2iMUiZZiRR;TO:Tl:S: 
Y ;  C;  77 ;  J71 ;  XI ;  X2  ;  DEL  %T  \A\T2\C1 
OBS+U . 0 
J+ 1 

T 1  -*-r  2  +  5)  (  N ,  M  )  p  2VZ1 0  <-Mp  0 
BO+BN+ 2p0 

X1+(M ,N)p (MxN)  NORRAND  1  0.5 
X2+(M,N)p(MxN)  NORRAND  2  1 
W<r(M,N)o(M*N)  EXPRAND  1 
Ol«-*((0.2xXl)+0.3xX2) 

Y+®(JlxW*  (*T2  ) 

DEL<-Y<,®OBS 
C1<-+/C<-+/DEL 
Y<-Yl®OBS 
IV+i*T 2 )xj 
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V 

VSTATLU1 

V  R+STAT  KiDF',ANSiII;JJtKKiTR;MS 
Cl]  II+JJ+KK+1 

[2] 

C3 ]  DF+  5  10  20  30  40  50 

[4]  TR+  10  20  50 

C5]  AMS«-(<(pW)x<prJ?)).K,4)p0 

[6]  MSE+ (C(pPF)x(p  TR ) ) , 7 ) p  0 

[7]  LOiANSLIIi  il^MS+K  THESIS (DFZJJ1 ,TRZKK1 ) 

C8]  MSEZIIi1  +  U+fMS)*K) ,  ((  +  /M5C;  1  2  4]*2)*K) 

[9]  II^JI+l 

[10]  +LO*\ (pTR)ZKK+KK+l 

[11]  (pDF)2JJ+JJ+KK+l 

[12]  fl+AMS 
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